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1. Introduction 

Hadronic matter deconfines at temperatures above 155 MeV, where quark and gluons consti¬ 
tute the relevant degrees of freedon. Nevertheless, the theory remains non-perturbative because 
of the strongly interacting infrared sector, and the equation of state (EOS) of this quark-gluon 
plasma requires non-perturbative analysis using methods such as lattice QCD. The EOS is needed 
in modeling the hydrodynamic evolution of the quark-gluon plasma probed in relativistic heavy-ion 
collisions and in understanding the cooling of the early universe. Phenomenologically, the most 
interesting region for lattice QCD analysis is between T siml40 MeV and 400 MeV, i.e., above 
the validity of hadron-resonance gas models and spanning the temperatures probed in relativistic 
heavy-ion collisions including the ‘transition’ region around 154(9) MeV [1]. 


2. Lattice Details and the Trace Anomaly 

We summarize the calculation of EOS [4] by the HotQCD collaboration using 2 -|- 1 flavors 
of fermions with the HISQ/tree action. In this calculation, we varied jS = 10/g^ in the interval 
5.9-7.825 along a line of constant physics (EoCP) defined by the ss meson mass tuned to M^s ~ 
695 MeV and the light quark mass fixed at mj = ms(20 {Mj^ k, 160 MeV), as used in our previous 
calculation [1]. To control discretization effects, the calculation was done at four values of the 
temporal size = 6, 8, 10, and 12. The spatial size was fixed at 4 times N^, i.e., Na/N^ = 4. The 
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Figure 1: The continuum estrapolation of (a) r\frf, (b) (c) riM^ and (d) wo/rj, along with their 

continuum values [2, 3]. The data have been corrected for the mass mistuning described later in the text. 
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Figure 2: (a) The scale, a/ri, normalized by the two-loop beta function, /(j3), showing a fit and a spline 
interpolation to the data and (b) the non-perturbative j3-function, Rp, compared to the two-loop value. 


lattice update was done with the RHMC algorithm with mass preconditioning [5]. All the results 
are reported with the lattice scale set by r\ = 0.3106(14)(8)(4) fm [6]. We checked that in our 
calculation this corresponded to ro = 0.4688(41) fm and wq = 0.1749(14) fm consistent with the 
previous determinations 0.48(1)(1) fm [7] and 0.1755(18)(4) fm [8] respectively. This estimate is 
also consistent, in the continuum limit, with the scale obtained from various fermionic quantities 
as shown in Fig. 1 . 

After integrating out the fermion degrees of freedom, the QCD partition function on an N^Nt: 
hypercubic lattice can be written as 

Zip,Nr) = f ^ ( 2 . 1 ) 

X.fi 


where sciU) and SfiU) are the gauge and fermionic actions in terms of the SU(3) link variables 
U. From this we can calculate the energy density, e, and pressure, p, in terms of the temperature, 
T, starting with the trace anomaly, 0^^ and using the relations 


^ dT \ T^ 


e — 'ip 
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p4 
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p4 

0^(7’) 

p4 




= -RpRm [2mi ((w)/,0 - ( W)/,t) +msi{\l/¥)s,o - (w)vr)] > 


( 2 . 2 ) 

(2.3) 

(2.4) 


where the subscripts T and 0 refer to expectation values at finite and zero temperature respectively, 
and I and s refer to the light and strange quark condensates, respectively. The function Rp is the 
non-perturbative j3 function, 

' (2.5) 


( 2 . 6 ) 


a \ dp 

shown in Fig. 2 and Rm is the mass renormalization function, 

1 dm,(j8) 


RmiP) = 


nisiP) dp 
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Figure 3: (a) The measured value of normalized by 695 MeV showing a few percent mistuning at higher 
j3. (b) The one-loop renormalization-group-invariant strange-quark mass in units of and the same after 
correcting for the mistuning, along with a spline interpolation and the asymptotic value at infinite j3. 


Since the measured value of has small deviations from LoCP, as shown in Fig. 3(a), we cor¬ 
rected it using lowest order chiral perturbation theory. The corrected mass parameter used in the 
analysis and to obtain R,„ is shown in Fig. 3(b). 

All our data for the trace anomaly are shown in Fig. 4(a) and compared to our previous cal¬ 
culations using the p4 and asqtad lattice actions. The noticeable difference is that the peak in 
the HISQ/tree data is lower and shifted to the left. In Fig. 4(b), we show that the results for 
r < 145 MeV at different Ni are consistent and agree with those obtained using the hadron reso¬ 
nance gas (HRG) model. 


3. Continuum Extrapolation 

The quality of the lattice data is sufficiently good that various strategies for extrapolating them 
to the continuum limit agree within errors. The one with least free parameters is a simultaneous 




(b) 


Figure 4: (a) The trace anomaly calculated at four different Nt compared to our previous calculations [9, 
10, 11]. (b) A magnification of the low temperature region with the results from a HRG model including all 
hadrons in the Particle Data Book [3] with masses less than 2.5 GeV shown as a solid line. 
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(a) (b) 

Figure 5: (a) A B-spline basis with two internal knots that are constrained to go to zero at 130 MeV. (b) 
The extrapolation of the lattice data at N-r = 8, 10 and 12 to the continuum is shown as a black band. The fit 
results at the various Nt are shown as colored bands. 


fit that interpolates the data in T and extrapolates in 1/A^^ —)• 0. This is implemented by choosing 
a basis of cubic splines that enforces the continuity of the function along with its first and second 
derivatives in T. An example of a B-spline basis with two knots is shown in Fig. 5(a). 

We explore fits with coefficients of the basis functions chosen as linear or quadratic functions 
of l/N^, thus enforcing the same continuity requirements on the extrapolation. To stabilize our 
fits, we match the value and the first temperature-derivative of the continuum extrapolated value 
on to the HRG results at T = 130 MeV. Since with n internal knots, the B-spline basis has n + A 
basis elements, the HRG condition at 130 MeV reduces the number of continuum coefficients to 
n + 2. For the extrapolation linear in 1 /V^, we have an additional n + A coefficients and the choice 
of n knot positions, for a total of 3n -|- 6 degrees of freedom (d.o.f). The fit was performed by 
minimizing the uncorrelated and n was chosen according to the Akaike Information Criterion 
(AlC) of minimizing 2 x d.o.f. All the calculations were done using the statistical program 
R [12] and its Hmisc package [13]. 

To estimate the errors in the fit parameters, we generated 20,001 synthetic data sets drawn 
assuming a normal distribution with variance given by the estimated error on each measurement. 
We assigned an independent conservative 10% error on both the value and slope of the HRG at 
r = 130 MeV. To the bootstrap error we linearly added a 2% error to account for scale uncertainty 
in determining T. 

An extrapolation ansatz just linear in 1 /N'^ does not fit the Nx = 6 data well. Our best estimate 
is obtained using a fit with 2 knots and linear in 1/V| to the = 8, 10 and 12 data at T < 400 MeV. 
To stabilize the fit at the upper end, we included the = 8 data up to T = 610 MeV. This final 
extrapolation is shown in Fig. 5(b). 
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Figure 6: (a) The continuum extrapolated integration measure (band) compared to the HRG estimate (solid 
line), (b) The continuum extrapolated pressure, p, energy density, e, and entropy density, s, (colored bands) 
compared to the HRG estimate (colored solid lines) and the Stephan-Boltzman value (solid black line) for 
an ideal gas. 

4. Results and Conclusions 

As shown in Fig. 6, our extrapolated continuum results for the integration measure, pressure, 
energy density and the entropy density agree in value with the HRG estimates for T < 150 MeV. 
The speed of sound and the specific heat shown in 7(b,c), which depend on higher derivatives of the 
partition function, are, however, not well predicted by this model. Figs. 6(b) and 7(c) also show that 
the quark-gluon plasma is still significantly far from a non-interacting gas even at F ss 400 MeV. 

In Fig. 7(a) and (b) we comapre our results with those by the Wuppertal-Budapest collabo¬ 
ration using a stout dicretization of the fermion action [14] and find good agreemenf. There may 
be a small discrepancy for T > 350 MeV which is probably unimporfanf for currenf heavy ion 
phenomenology, buf needs furfher invesfigafion fo sfudy fhe approach fo fhe high temperafure per- 
furbafive regime. Lasfly, fhe disagreemenfs wifh our previous results for the EOS using the p4 and 
asqtad actions were due to large cutoff effects in those calculations. 
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Figure 7: (a) Comparison of our results for pressure, p, energy density, e, and entropy density, s, (colored 
bands) with the results from the Wuppertal-Budapest collaboration using the stout discretization [14] (grey 
bands), (b) Our results for the speed of sound, Cj (colored band) compared to the stout results (data points) 
and the HRG estimate (black solid line), (c) Our results for the specific heat, Cy, (colored band) compared 
to the HRG results (colored line) and the value for a non-interacting gas (horizontal black line). Also shown 
is the result from a parameterization designed to the fit the pressure and energy density (black curve). 
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